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Abstract 

A model independent analysis has been performed on the LEP and SLD data 



on Z decays. Using only very weak theoretical assumptions, the effective vector 
and axial-vector couplings of leptons, c quarks and b quarks have been extracted. 
Although the lepton and c quark couplings agree well with Standard Model predic- 
tions, those of the b quark show deviations of more than three standard deviations. 
The effect is mainly in the right-handed b quark coupling, the left-handed cou- 
pling being consistent (at the 2a level) with the Standard Model prediction. The 
probability that the observed deviations of all the measured effective couplings are 
statistical fluctuations from lepton universality and the Standard Model is estimated 
to be 0.9%. The estimated probability that the deviations in the leptonic and b 
quark couplings alone are a fluctuation is 0.18%. A thorough discussion is made of 
the internal consistency of the different measurements contributing to the average 
values A[ and used to extract the b quark couplings, as well as possible sources 
of systematic error that may not, hitherto, have been taken into account. Excluding 
r-polarisation measurements, which show internal inconsistencies, from the aver- 
ages increases the deviations of the extracted b quark couplings from the Standard 
Model predictions to the four standard deviation level. 
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1 Introduction 



This study is based on a recent compilation |]J of experimental results on Z decays. 
The aim is to answer the question: 'are the data consistent with the predictions of the 
Standard Electroweak Model (SM) [fj]' To this end the analysis is carried out in three 
independent steps. In the first, the data is used to extract the effective vector (v) and axial 
vector (a) coupling constants of the charged leptons (assuming lepton universality) and of 
the c and b quarks. This is done using only weak theoretical assumptions. The effective 
couplings are then compared with the predictions of the SM and the confidence levels 
(CLs) for consistency of the measurements with the SM are calculated. In the second 
step the experimental data contributing to deviations observed from the SM in the b 
quark couplings are critically examined. Issues addressed include the roles of statistical 
and systematic errors, as well as the internal consistency of physical parameters measured 
using different experimental methods. Finally the observed deviations from the SM are 
assumed to represent real physical effects, whose interpretation is discussed. The three 
steps described above constitute the material of the following three Sections. A summary 
and outlook are given in the final Section. 



2 Extraction of the Effective Weak Coupling Con- 
stants 

It is convenient to define the following auxiliary quantities^: 

rf = v f /af (2.1) 
*/ = («/) 2 + (^/) 2 (2-2) 

that may be simply derived from the measurements. The experimental errors on ff and Sf 
(f here stands for lepton or quark) are, unlike those in Uf and a/, essentially uncorrelated, 
simplifying the calculation of the statistical significance of any deviations observed from 
the SM expectations. Throughout this Section, the SM predictions quoted are those of 
the global SM fit with m t = 172 GeV, m H = 149 GeV reported in Ref.[l]. The effect 
of varying the Higgs boson mass, the only remaining unknown parameter of the SM, is 
discussed in Section 4 below. 

The quantities ff (/ = I, c, 6)0 and s; may be directly obtained from the data without 
any additional assumptions concerning the poorly-measured couplings of the u, d, s 
quarks. A further assumption is, however, necessary in order to extract s c , and hence 
the c and b quark couplings. In order to perform an analysis which is, as far as possible, 
'model independent', and to avoid the specific assumption of the validity of the SM, as 
used in the fits of Ref. [1], the weaker hypothesis of quark- lepton universality is made for 
the fermions e, /i, r, u, d, s, c. That is, all these fermions are assumed to have the same 
effective weak mixing angle []]]. The derived values of v c and a c presented below are found 

1 The fermion masses are set to zero in Eqn.(2.2). Only for the b quark do the fermion mass terms 
give a non-negligible contribution 

2 Unless otherwise stated, I is a generic lepton label and e — /i — r universality is assumed 
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to be, within errors, in good agreement with this hypothesis. Another possibility is to 
assume a value of a s (Mz) derived from non electroweak- related measurements. In this 
case the c and b quark couplings may be derived from the ratios Tq/Ti ( Q = c,b) without 
any assumption concerning the couplings of the light quarks. There is now, however, the 
disadvantage that the extracted values of the electroweak couplings are strongly correlated 
with the assumed a s (Mz) value. 

At LEP, ff is found from the measured, corrected, pole forward/backward charge 
asymmetries Ap B via the relations: 

A% = \A e A f (2.3) 

f - (a f y + (v f y i + r) [ ■ } 

A e and A T have also been measured at LEP via the angular dependence of the r- 
polarisation asymmetry: 

where 

F{9) = 2 cos 0/(1 + cos 2 6) 

and 9 is the angle between the incoming e~ and the outgoing r~ in the r-pair centre 
of mass frame. At SLD, A e is directly measured by the left/right beam polarisation 
asymmetry Alb.-, while A c and Af, are determined from the left /right-forward/backward 
asymmetries of tagged heavy quarks. 

The separate LEP and SLD average values of the electroweak observables, which 
are directly sensitive to the effective couplings, are reported in Table 1. The combined 
LEP/SLD averages of A h A c , A b , R c and Rb, where Rq = YQ/Y had , ( Q = c, b) are 
reported in Table 2. 

It may be remarked that, while there is good agreement between the different values of 
Ai derived from A°/ B (I = e, /i,r) using Eqn.(2.3) and that derived from A LR , ( weighted 
average 0.1533(27)0, x 2 = 3.85 for 3 DOF , CL = 28%), the values of both A e and A T 
derived from the r-polarisation measurements are significantly lower. In fact, the average 
value of A\ from r-polarisation lies 2.5 a below the average from A°p l B , A^ B . Including, 
or not including, the r-polarisation data changes the LEP/SLD average value of Ai by 
more than one standard deviation (see the first column of Table 2). Because of this 
possible inconsistency in the measured A\ values from different sources, the extraction of 
the coupling constants will be done throughout this paper using values of Ai that either 
include, or exclude, the r polarisation results. Significant differences are found. Possible 
explanations for the apparent inconsistencies in the Ai measurements are discussed in the 
following Section. 

The values of ff (f = l,c,b) derived from the measured values of A\ using Eqn.(2.4) 
are presented in Table 3. For the b quark, mass effects were taken into account by using 



3 Throughout this paper total experimental errors are given in terms of the last significant figures. 
0.1533(27) denotes 0.1533± 0.0027 
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the corrected form of Eqn.(2.4): 



2(yi - 4/i b )r fe 
1 - 4/x 6 + (1 + 2/i fe )r^ 



M = , T n (2-6) 



where fi b = (m b (M z )/M z ) 2 ~ 1.0 xlO -3 . The running b quark mass is taken as rrib(Mz) = 
3.0 GeV 0. Agreement is seen with the SM at the 2a level for Fj, at < lcr for F c , but only 
at the 3.3a level for F&. The similar discrepancy for A b was mentioned, but not discussed 
in terms of the b quark couplings, in Ref. [1]. If the r-polarisation measurements of Ai are 
excluded from the average, the discrepancy of r\ with the SM approaches 3aQ and that of 
F& exceeds 4a. 

The quantity si is derived from the leptonic width Ti using the relation: 

_ /_ \2 /— \2 127lT2 1 . . 

* = {T " ] + w = v^M (i + Mg^) (27) 

The value obtained for si, quoted in the first column of Table 4, uses the LEP average 
value of from Table 1 together with: G M = 1.16639 x 10~ 5 (GeV) 2 [|, M z = 91.1863 
GeV, and a(Mz)^ 1 = 128.896 JT|. Good agreement is found with the SM value. Solving 
Eqns.(2.1) and (2.2) for ai and vi yields the results presented in Table 5. As in the 
calculation of all the other effective couplings, the signs of Hi and vi are chosen to be 
the same as the SM predictions. The values of Hi and vi are in good agreement with the 
LEP+SLD averages quoted in Ref.[l], taking into account the slightly different analysis 
procedures^. Both ai and vi are in agreement with the SM predictions. 

The quantities sq ( Q = c,b), including quark mass effects, may be derived from the 
measured quantities Rq via the relation: 

SQ = («Q) 2 (1 - 6/ig) + {v Q f = - - R fo Q ED„QCD (2-8) 



;i - R Q )c^c t 



where 



and 



Sq= E[(s) 2 ( 1 -6/i,) + 



v 



C% = l + 5 i Q -<5 i g ^ Q > (i = QED,QCD), N = for q ± b 
6 qed = !h£ a / Mz)) 5 QCD = L0Qas + 1A2a 2 s? OD = .99a s -lMa 2 s 

4:71 ^ 

q is a generic quark flavour index, e q the quark electric charge in units of that of the 
positron and a s = a s (Mz)j^ ■ < X > denotes the quark flavour average of X. As 
mentioned above, fib = 1.0 X 10~ 3 while, taking into account the present experimental 
error on R c , fi c is set to zero. The numerical values of the QED and QCD correction 

4 As shown in Section 4 below, the discrepancies with the SM predictions for the leptonic couplings, 
unlike those of the b quark, are reduced by assuming a smaller value of m# and a larger value of m t than 
that found in the global fit of Ref. [1] 

5 Ref. [1] included small mass corrections in calculating a; and vi which are neglected here. 
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factors, with a s (Mz) = 0.12 and a(Mz) 1 = 128.9, are presented in Table 6. The non-b 
quark couplings in Eqn.(2.8) are written, conventionally, as: 

a q = Jp- q Tl (2.9) 
v q = v^( T 3 - 2e 9 (^) 2 ) (2.10) 



where, assuming non-b quark lepton universality []: 

yfc = y/pi = 2\a t \ (2.11) 
(si)" = 1(1 -n) (all q^b) (2.12) 

and T3 is the third component of the weak isospin of the quark q. Substituting the 
measured values of ri, ai, from Tables 3, 5 and of R c , Rb from Table 2, leads to the 
values of s c , ~s~b reported in Table 4. Note that the value of ~s~b, and hence Tib and Ub are 
extracted first. The latter are then substituted into Eqn.(2.8) (taking into account their 
experimental errors) in order to find s c . In Table 4 good agreement is seen between the 
measured values of ~s~i and s c and the SM predictions. On the other hand, ~s~b lies 1.3 a 
above the prediction, a residual of the well known L Rb problem' [I]]. Solving Eqns.(2.1) 
and (2.8) then gives the effective coupling constants for the heavy quarks presented in 
Table 7. The values found, as well as the errors, agree well with those reported by Renton 
in a recent review |7j. The solutions for af,Vf obtained from the essentially uncorrelated 
quantities Tf and Sf are shown graphically in Figsla,lb,lc for f = /,c,b respectively. The 
corresponding solutions when the r-polarisation measurements of Ai are excluded from 
the average are shown in Figs.2a,2b,2c. It is clear from Figslc,2c that largest discrepancy 
with the SM is in the parameter ^(completely determined by Ab) rather than in Sb 
( essentially determined by R b ). Indeed, if the SM value for the latter is used, instead of 
the measured one, to solve for a b and Vb, the discrepancies between the values found and 
the SM are almost unchanged. 

Although the c quark couplings agree well with the SM, and are also consistent with 
the quark-lepton universality hypothesis, both Ub and Vb differ from the SM values by 
more than three standard deviations. The errors on these quantities are, however, highly 
correlated. The statistical significance of these deviations is discussed in detail below. 

It should be remarked that, although a particular value (0.12) of a s (Mz) has been 
assumed in order to extract the effective couplings of the heavy quarks, the sensitivity 
to the chosen value is very weak. Varying a s [Mz) over the range 0.1 < a s [Mz) < 0.14 
leads variations of only ~ 3 x 10~ 4 in Ub and Vb to be compared with experimental errors 
~ 1 -4 x 10~ 2 (see Table 7). 

A further constraint on the quark couplings is provided by the measurement of the 
mean quark forward/backward charge asymmetry: 

{ FB) ~Y. q [{i-^ q )(a q y + {v q y] (2 ' 13) 



6 Here the weak isospin symmetry of the SM is invoked to calculate the unobserved couplings. It is 
also assumed that the quantum corrections contained in p q and (si) 2 , though not necessarily those of 
the SM, are universal. 
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Figure 1: Constraints on the effective couplings a/, Vf provided by the measurements 
of Tf and Sf . a) leptons, b) c quarks, c) b quarks. The cross-hatched areas show ±ler 
limits. The dotted lines in a),[c)] show 2a, [3a] limits for r;,[r&]. SM is the Standard 
Model prediction for m t = 172 GeV, m H = 149 GeV. 
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Figure 2: As Fig 1, except that r-polarisation measurements are excluded from the LEP 
average value of A\. The dotted line in a) shows the 2a limit for r 7 /. The dash-dotted 
[dotted] line in c) shows the 3a, [4a] limit for r&. 

All experimental analyses performed to date have assumed the correctness of the SM 
and have used measurements of (A q FB ) to determine a value of sin 2 O^fj [fl- Inserting 
the average value of the latter reported in Ref. [1] into the SM formula for (A q FB ) and 
propagating the error leads to the 'measured' value: 

(A FB ) = 0.1592(86) 

As shown in Table 8 this value is consistent with the SM prediction, with the 'model 
independent' prediction given by inserting the lepton and b quark couplings from Tables 
5 and 7 into Eqn(2.13) and assuming non-b quark lepton universality for the u,d,s,c 
quarks, as well as the prediction when, in the latter case, the measured b quark couplings 
are replaced by the SM ones. With the present experimental errors, (A q FB ) is therefore 
insensitive to possible deviations of the b quark couplings from the SM, of the magnitude 
observed in the A^ measurements. 

As mentioned earlier, in order to avoid having to introduce an accurate value of a s (Mz) 
as a correlated parameter in the extraction of the heavy quark effective couplings, the 
hypothesis of non-b quark lepton universality was made in deriving the value of s& from the 
measured quantity R^. The consistency of this assumption may be checked by extracting 
a s {Mz) from the LEP average value of R\ = Thad/^i 

Ri = 20.778(29) 
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using the relation: 

<C QED ><C QCD >E 

Ri = 3 q - 2.14 

c qed 5j v ; 

The QED and QCD correction factors < > and < C® > are averaged over all 

quark flavours. The QED correction factors are: 

< C® ED >= 1.00040, Cf ED = 1.0019 

Inserting the measured values of Sf, and s;,and using non-b quark lepton universality to 
evaluate s q ( q ^ b ) gives, for the QCD correction factor: 

< C® CD >= 1.0394(21) 

Using the third order perturbative QCD formula |§: 

< >= 1 + im as{Mz) + 0.9(^M) 2 - 15(^M) 3 (2.15) 



gives: 

a s (M z ) = 0.116lffi 
which may be compared to the global fit value of Ref. [1]: 



a s {M z ) = 0.120(3) 

The good agreement of the model independent analysis result with the global world av- 
erage value: a s {Mz) = 0.118(5) found in two recent reviews ||, [L0| of all published 
measurements of a s , shows that an analysis assuming this value of a s (Mz), but without 
the assumption of non-b quark lepton universality, would lead to essentially the same 
values of the b quark couplings as those reported in Table 7. In the fit used in Ref. [7] 
to determine the heavy quark effective couplings the constraint a s [Mz) = 0.123(6) was 
imposed. As mentioned above, the fitted heavy quark couplings are very consistent with 
those found in the present analysis. 

In order to correctly calculate the statistical significance of the deviations from the SM 
predictions of the effective couplings shown in Tables 5 and 7 it is necessary to take into 
account the correlations between the errors of the different quantities. To avoid the very 
large correlations between the errors on a/ and Vf (for the case of b quarks the correlation 
coefficient is -0.96) it is convenient to use, in calculating the x 2 , the equivalent quantities 
7f, Sf for which the errors are uncorrelated for a given fermion flavour /. Important 
correlations still exist, however, between the errors on (rj, r c ) and (Fj, r&) in the case that 
r c and Fj are extracted from forward/backward asymmetries using Eqns.(2.3),(2.4) and 
(2.6). The correlation coefficient is: 

ClQ = Jlz^±^l^X (Q = C , 6 ) (2 . 16) 

Substituting the parameters from Table 3 gives: 

C lc = -0.29, C lb = -0.52 
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The results on the CLs for the agreement with the SM of different sets of effective weak 
coupling constants, parameterised in terms of ff and Sf, are collected in Table 9. These 
CLs assume perfect statistical consistency of the different measurements contributing to 
the averages. The entries in the first column of Table 9, giving the level of agreement of 
(fi, si) with the SM prediction are simply calculated from the entries of Tables 3 and 4 
using a diagonal error matrix, since the errors on fi and s; are uncorrelated. Calculating 
separately the contributions to x 2 from r\ and r&, where the latter is derived from the LEP 
A°p B measurement, and Tb derived via Eqn.(2.6) directly from the SLD Aj, measurement, 
gives the entries reported in the second column of Table 9. The CL for agreement with the 
SM prediction of 1.4% drops to only 0.06% if the T-polarisation measurements of Ai are 
excluded. The third column of Table 9 results from adding to the x 2 i n the second column 
the (uncorrelated) [] contributions of s/ and Sb- In the fourth column of Table 9 the x 2 an d 
CL of the variables Ti, Tb and T c taking into account the Ti-Tb and T\-T c correlations are 
given. In the last column of Table 9 the (uncorrelated) variables si, Sb and s c are added 
to those of the fourth column. Note that the number of degrees of freedom corresponding 
to the x 2 values reported in the second, third, fourth and fifth columns of Table 9 are 
3, 5, 5 and 8 respectively, since the T c and T b measurements derived from the SLD A c , 
A b determinations give separate, uncorrelated, contributions to the x 2 ■ As expected, the 
agreement with the SM improves as the number of degrees of freedom of the x 2 increases 
(the more parameters are considered, the more likely is a deviation associated with any one 
of the parameters to be consistent with a statistical fluctuation). However, there is still a 
factor ~ 10 difference between the CLs including (or excluding) the r-polarisation data. 
Taking into account the CL (8.4%) for self- consistency of the different Ai measurements, 
the probability^ that all six effective couplings are consistent with lepton universality and 
the SM is 0.9%. The similar probability for the leptonic and b quark couplings alone is 
0.18%. If the r-polarisation measurements are excluded the latter probability drops to 

0. 018%. 

It is important emphasise that a correct calculation of x 2 and the associated confidence 
levels requires that all relevant correlations between errors are taken into account. If a 
X 2 is calculated from 'raw' experimental measurements, such as those presented in Table 

1, erroneous conclusions as to the consistency of the data with the SM will be drawn. 
Assuming non-b quark lepton universality, the 14 measured electroweak observables pre- 
sented in Table 1 depend on only four unknown parameters, the effective couplings of the 
leptons and of the b quarks. The 'raw' x 2 calculated from the 'pulls' |1J in the last column 
of Table 1 is 21.3 for 14 dof (CL = 0.093 ) or, excluding the r-polarisation data 19.5 for 
12 doff\ (CL = 0.077 ). The strong sensitivity of the CL to inclusion or exclusion of the 
r-polarisation data is completely lost using the 'raw' x 2 ■ m f ac ^ the ~ 2er effects seen in 
the 'pulls' of the observables A°p T B , A e (SLD), A°p b B and Ab (SLD) add, because of corre- 
lations, constructively in the parameter Tb to give the observed deviation from the SM of 

7 Actually there is a weak correlation between Sb and r/ following from Eqn.(2.8), where r; is used to 
calculate Sq. However the correlation coefficient is only ~ 0.08 and is neglected here. 

8 Here the term 'probability' is used in the the usual sense of the fraction of all cases expected to 
have a CL less than the observed value. For independent \ 2 tests the probablities are assumed to be 
uncorrelated 

9 In Rcf [1] , the SM prediction is obtained by fitting m t , m# and several other electroweak parameters 
to the observables of Table 1 as well as others which are not directly sensitive to the effective couplings 
(see Table 20 of Ref.[l]). Here, for comparison purposes, the values m t — 172 GeV, m# = 149 GeV are 
assumed so that the SM prediction has no free parameters 
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> 3a. First extracting the essential theoretical parameters (the effective couplings), and 
then comparing with the SM predictions using a x 2 test, is more sensitive to deviations of 
these parameters from the SM predictions than the 'raw' x 2 ■ For the latter an inevitable 
statistical dilution occurs because of the large number of experimental observables (14) 
used as compared to only four effective coupling constants that determine the theoretical 
prediction. This dilution effect becomes even more marked when additional observables, 
not directly sensitive to the effective couplings, are added to the x 2 as i n the global fit 
of Ref.[l]. In fact, the five additional observables used in the fit contribute only 0.55 to 
the total x 2 °f 19.1, indicating an overestimation of the errors on these quantities^ that 
reinforces the statistical dilution. 



3 Discussion of the Experimental Measurements of 
Electroweak Observables 



In order to derive the average values of the electroweak observables presented in Tables 
1 and 2, many different experimental measurements were combined [ffl]. In the light of 
the apparent deviations seen from the SM predictions in the model independent analysis 
described in the previous Section, the first question that should be asked is whether the 
experimental measurements are reliable and consistent. An important general question 
is whether the uncertainty on the measured observable is dominated by statistical or 
systematic errors. Only in the former case, the error can be interpreted, with confidence, 
in the statistical sense (la = 68% CL) and the probabilistic meaning of the CL of a x 2 
test can be expected to be reliable. This is no longer the case if the systematic error 
is dominant. As there is no definite, agreed, procedure for assigning systematic errors, 
the meaning of the error can depend on psychological (or even sociological) factors. If 
the physicist is over-conservative in assigning the error, real deviations from a theoretical 
expectation can be missed, or in the contrary case spurious detector related effects wrongly 
interpreted as 'new physics'. Some check on the degree of conservatism, or otherwise, of 
physicists is however provided if there are many repeated measurements of the same 
quantity, as provided, in the present case, by the different experiments at LEP. When the 
errors are dominated by systematic effects, or contain a large systematic contribution, 
then comparing the weighted average error to that calculated by applying the Central 
Limit Theorem to the different measurements of the same quantity gives an indication 
whether the systematic errors are over- or under-estimated. Such a test only applies to 
errors which are uncorrelated for the different measurements of the same quantity. This 
test may be best applied when each experiment measures the same observable using a 
similar method. If consistency is found in this case, but inconsistencies are found when 
the same physical quantity is measured using different observables ( for example, the 
quantity A\ may be measured using either forward/backward asymmetries for different 
dilepton final states, or r-polarisation ) it is probable that there is an unknown source 
of correlated systematic error. As discussed further below, such an error can arise, for 
example, due to an inadequate treatment of QED radiative corrections. 

The apparent deviations from the SM predictions of the b quark couplings seen in 
10 A x 2 °f 0.55 for 5 degrees of freedom correponds to a CL of 0.990 
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Table 7 are due those observed in only two of the basic electroweak parameters: At, and 
A\ . Most of the current information on At, is derived from the Z-peak A°p b B measurements 
of the four LEP experiments. The error on this At, measurement is roughly half that of the 
SLD At, measurement. The individual measurements of the experiments [Q contributing 
to the LEP average A°p b B measurement quoted in Table 1 are shown in Fig. 3, together 
with the weighted average value, its error, and the SM prediction. There is no hint of 
any badly understood systematic effect in the distribution of these measurements. Except 



for the ALEPH jet-charge measurement [EL], all errors, even those using all LEP1 data 



are statistically dominated. The weighted average value agrees well with the individual 
measurements (x 2 = 5.9 for 8 dof, CL= 66%). The estimate of the error on the mean 
value, given by the sample variance of the seven most accurate determinations, using 
the Central Limit Theorem is 0.00232, in excellent agreement with the weighted average 
error of the same data points, which is 0.00237. The situation is very different for the 
parameter A\. Shown in Fig. 4 are the values of A\ derived from the LEP measurememts 
of A°p Bl A°p B , A°p T Bl A e and A T from r-polarisation and A e as measured by SLD using 
the left/right beam polarisation asymmetry. Although the overall consistency of the 
individual measurements with the weighted average value seems acceptable (x 2 = 9.7 for 
5 dof , CL= 8.4%) the internal consistency of the various measurements is much worse. 
In particular there are three, essentially independent, Q ~ 2 — 3a deviations concerning 
r-related measurements: 

• A°p T B is 1.8cr higher than the average of A°p B and A°p B . Also, in all 4 LEP exper- 
iments, (see Table 3 of Ref.[l]) A°p B is higher than A°p B or A°p B . Assuming no 
systematic bias, the probability for this is 1 in 81. 

• The average value of Ai extracted from the r-polarisation data: 0.1393(50) lies 2.2<r 
below that, 0.1522(30), given by the weighted average of the non-r measurements. 

• Assuming lepton universality, the LEP average value of Ap B gives, using Eqn.(2.3), 
A\ = 0.1649(71). This is 2.9cr higher than the mean A\ calculated from the r- 
polarisation measurements of A e and A T . 



These deviations cannot be explained by a breakdown of lepton universality for the r. 
Using the r-polarisation A e measurement to extract, using Eqn.(2.3), A T from A% gives 
A T = 0.1968(204). This may be compared with the r-polarisation measurement: A T = 
0.1401(67). There is a 2.6<r discrepancy. In the case of a breakdown of lepton univer- 
sality, the two determinations of A T must give a consistent result that is significantly 
different from the measured A e value. In fact, the values of A e and A T found using 
r-polarisation are consistent within 0.19<r. 

Assuming gaussian errors, the probability that all these r-related apparent deviations 
from lepton universality are statistical fluctuations, is 7.5 x 10~ 6 . This situation may 
be compared with that for the non r-related measurements of A\ derived from A°p B , 
Ap B and A LR which are, respectively 0.147(11), 0.148(61) and 0.154(4). Labelling these 

11 There is a weak correlation between two of the three effects ,in that the A°p C B , A° FB measurements, 
used in the first consistency check, contribute also to the 'non-r' weighted average value of Ai used in the 
second. This can be avoided by comparing the r-polarisation value of Ai with the Alr measurement. In 
this case an even larger discrepancy of 2. 4a is found. 
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Figure 3: LEP on-peak A°p b B measurements. ALEPH: solid squares, DELPHI: open 
circles, L3: open triangle, OPAL: open squares. The ±ler region around the weighted 
average value is indicated by the hatched band. The vertical line is the Standard Model 
prediction for m t — 172 GeV, ran = 149 GeV. 
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Figure 4: LEP and SLD Ai measurements. The hatched band shows the ±1<t re- 
gion around the weighted average value. The weighted average value, excluding the r- 
polarisation measurements, is given by the dashed vertical line. The solid vertical line 
is the Standard Model prediction for m t = 172 GeV, itlh = 149 GeV. The open square 
shows the value of Ai derived from the LEP average value of A°p b B assuming the SM; this 
datum is not included in the weighted averages shown. 
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measurements 1,2,3 respectively, the deviations between the pairs 1-2, 1-3 and 2-3 are 
0.08a, 0.6a and 0.83a. The measurements are perfectly consistent. 

It may be remarked that the A°p B measurement, which occurs in the first of the three 
above-mentioned apparent deviations, has been included in the evaluation of the average 
values of A\ used in the model-independent analysis of the previous Section. However, 
excluding this datum from the Ai average in addition to the r-polarisation data gives 
A\ = 0.1516(29). Thus, (see Table 2) the deviations from the SM predictions will lie 
between the 'r-polarisation in' and 'r-polarisation out' cases discussed above, if all 
r-related measurements are excluded. 

Also shown in Fig. 4 is the value of Ai derived from A°p B , assuming the correctness 
of the SM. The value so obtained, 0.1396(33), differs from the weighted average of the 
purely leptonic measurements by 2.6<r, or by 3.2cr if the r-polarisation data are excluded. 
It is clear, from the analysis of the previous Section, that these discrepancies are mainly 
due to the deviations of the b quark effective couplings from the SM predictions. The 
quantity sin 2 9 l ^j used in Ref. [1] is directly related to A\ via Eqns.(2.4),(2.12). The poor 
consistency of the different sin 2 9 l ^jj determinations in Table 19 of Ref.[l] is largely due 
to the inclusion of values derived from A°p B and (A q FB ) assuming the correctness of the 
SM. The common origin, in the b quark couplings, of the poor agreement of the different 
sin 2 determinations and the 3<r deviation of the measured LEP-SLD average value of 
At, from the SM prediction, was not pointed out in Ref.fl]. 

Each type of observable contributing to the average values of Af. A°p l B , (I = e,/i, r) 
A e and A T from r-polarisation and A e from Alr are now discussed it turn, taking into 
account the internal consistency of measurements of the same quantity performed by 
different experiments, the relative importance of statistical and the estimated systematic 
errors, and possible systematic effects that may not have, so far, been taken into account. 

The LEP Leptonic Forward/ Backward Charge Asymmetry Measurements 

The values of A°/ B , (I = e, fi, r) for the different LEP experiments are presented in Table 
3 of Ref. [1]. A detailed breakdown of the systematic errors of the different data sets is 
found in Table 2 of the same paper. Although the individual experimental arrors are 
usually statistics dominated, the statistical and systematic contributions to the error on 
the weighted average (see Table 20 PI of Refill) are almost equal. The most remark- 

A I 

able systematic feature of the A FB measurements, as already mentioned above, is the 
relatively high value of A FB found by all four experiments. Since each of the LEP exper- 
iments have comparable statistics in each channel, each can be considered to provide an 
independent estimate, with a similar weight, of A°/ B . Disregarding the estimated errors 
on each measurement, an unbiased statistical test of the consistency of the measured val- 
ues of A°p B can be made using the Student's t distribution. This is done by calculating 
the probability that, say, the A°p T B and A°p B measurements of the different experiments 
are consistent with a common mean value |12|| . The results of this comparison for the 3 
possible pairings (r, /i), (r, e) and (/x, e) are presented in Table 10. Also shown in this 
Table are the CLs for consistency of the measurements, based on the total errors quoted 
in Table 1. The systematically larger CLs found using these errors is perhaps indicative 

12 Note there is a mis-print in this Table. The estimated systematic error on A°p B should presumably 
be 0.0007 not 0.007. 
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that the point-to-point systematic errors tend to be over-estimated. The good agreement 
between the /i-e measurements and the poor agreement between t-\i and r-e indicates the 
possible presence of a correlated systematic effect, not included in the present systematic 
error estimate, for the r-pair channel. 

An obvious candidate for for such an effect is the QED radiative correction. On the 
Z-peak this is large; the combined 0(a) and 0(a 2 ) corrections amount to ~ —110% of 
the corrected pole asymmetry for muon-pairs [|13|]. The systematic error on the weighted 
average value of A°/ B is then ~ 4.4% of the radiative correction, and the observed deviation 
of the r from the combined /i-e results is ~ 27% of it. It seems, however, unlikely that 
the error on the theoretical estimate of the QED radiative correction, essentially due to 
unobserved initial state radiation, and the associated virtual corrections, could be large 
enough to explain the high value of the r-pair asymmetry. One effect that can produce 
large changes in the forward/backward asymmetry is intial/fmal state interference in 
the case that hard cuts are applied to the radiated photons |L4], jlS] . Indeed, for some 
decay channels, such as ttu, -kh-kv, tight cuts are applied in some of the LEP analyses 
in order to cleanly separate them from channels containing extra 7r°s. Such effects could 
be investigated by comparing the forward/backward asymmetries for different r decay 
modes and different cuts on additional photons. 

The LEP t -Polarisation Measurements 
Results of the measurements of the average r-polarisation of the four LEP experi- 



ments [16, 17, 18, 19], for each r decay channel analysed, are presented in Table 11. 
The first error shown on each measurement is statistical, the second one systematic. The 
weighted average values A T and errors are given separately for each experiment and for 
each decay channel. The overall weighted average and its error are also given. The results 
shown in Table 11 do not correspond exactly to those used in the LEP averages reported 
in Ref.[l], although there is a large overlap. The data sets chosen are those for which 
both statistical and systematic errors have been given for each decay channel. 

It can be seen that, unlike for the forward/backward charge asymmetry measurements, 
the quoted statistical and systematic errors of individual experiments are comparable for 
almost every decay channel. The systematic error estimates can be tested by comparing 
the weighted errors, for each experiment or each decay channel, with the error estimate on 



the mean value given by the Central Limit Theorem: yJ2( x — x) 2 /n(n — 1). The latter 
error estimates are given in Table 11 in square brackets next to the weighted averages. 
It can be seen that the level of agreement is better between experiments than between 
decay channels. For the euu, fiuu, ttu, and a\v channels there is an indication of an over- 
estimation of the systematic errors, while for pv the error calculated from the sample 
variance is much larger than the weighted average error. This is due to the wide spread 
of the measurements in this channel; the DELPHI and L3 values are much larger than 
those of ALEPH and OPAL Q. 

The common mean test using the Student's t distribution, applied above to the for- 
ward/backward charge asymmetry measurements has also been applied to the measurents 
of A T presented in Table 11. The results of the comparisons of different decay channels are 



13 The weighted average of DELPHI and L3 is A T = 0.171(19), that of ALEPH and OPAL A T = 
0.113(15). The difference is 2.5 standard deviations 



14 



presented in Table 12 and of different experiments in Table 13. In each case the CLs for 
consistency based on the total experimental errors are also given. The main systematic 
features of the A T measurements are the relatively high values found by all experiments 
for the evv channel and the high value of the L3 weighted average, as compared to those 
of the other three experiments. The latter effect is also seen in the r-polarisation mea- 
surements of A e I] . The systematically larger value for the evv channel is reflected in the 
poor confidence levels (0.35%-2.5%) of the Student's t tests for the channels e — /x, e — n 
and e — ai in Table 12. Relatively worse CLs are also seen for these channels in the test 
based on the total errors, but as expected from the larger experimentally assigned errors 
as compared to those calculated from the sample variance, higher absolute CLs (11%- 
31%) are found. The agreement between the two types of test is much better in Table 
13, from which one might be tempted to conclude that all 4 experiments give consistent 
result s|3- 

As in the case of the r forward/backward charge asymmetry, the radiative correction 
associated with final state radiation appears as an obvious candidate to explain both the 
systematic differences observed between different decay channels and those, mentioned 
above, between the r-polarisation and other measurements of A\. The final state radiative 
corrections in the r-polarisation measurements are not only large for most decay channels, 
but depend strongly on the detection efficiency of the radiated photons and hence on the 
acceptance and resolution of the LEP detectors as well as the experimental cuts. The 
r-polarisation is measured by fitting the energy spectra of r decay products. The latter 
are directly effected by the rate and energy spectra of, and cuts applied to, the radiated 
photons. Any systematic errors in the treatment of final state radiation are thus directly 
correlated to systematic errors in the r-polarisation measurement. Consider, for example, 
the measurement of (P T ) using the nv decay channel on the Z peak. Measuring only the 
pion energy and neglecting that of the radiated photons has been estimated [Q to shift 
(P T ) by ~ 20% of its value. This may be compared to the systematic error assigned to 
the LEP average value of (P T ) of 3.3% of its value |IJ. Since radiative corrections are not 
included in the systematic error estimate, a tacit assumption is thus made that they are 
known to much better than ~ 10% of their value. 



All of the LEP experiments have used the same Monte Carlo program, KORALZ |2T 
to correct for radiative effects in fitting the energy spectra of the r decay products. No 
experimental checks (measurements of the rate and distributions of final state photons) 
have been published, and any systematic errors assigned by the experiments for radiative 
corrections have been very small in comparison to detector-related sources of systematic 
error. As for the case of the r forward/backward asymmetry it seems unlikely, however, 
that the approximations made in simulating the radiative corrections can account for the 
all the different values of A T found using different decay channels. For the decay mode of 
largest statistical power, r — > pv, where the sensitivity to the r-polarisation is given by 
angular information from the p decay as well as the energy distributions of decay products, 
the sensitivity to radiative corrections is small, The entire radiative correction generated 



14 The Student's t test is less well adapted to the comparison of the different LEP experiments since the 
different decay channels have very different statistical sensitivites, whereas each test datum is treated, 
in the test, on an equal footing. On the other hand, since each experiment has comparable statistics 
and sensitivity for a given decay channel, the Student's t common mean test is well adapted to the 
channel-by-channel comparisons 
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by KORALZ has been shown to generate a shift of only -0.011 in A T in this case0. For 
the case of the itv decay channel, a large 'structure dependent' effect in the final state 
photon spectrum due to the decay a\ — > tc-/ is to be expected. This contribution and 



its interference with the bremsstrahlung amplitudes, have been estimated [22] to change 
the yield of high energy photons with x = 2E 1 /m T > 0.6 by a factor > 2. However, 
the radiative correction to A T is dominated by soft bremsstrahlung photons, so no large 
corrections are to be expected from this effect. For the euu, [ivv decays KORALZ uses 
the exact 0(a) matrix element. Particularly for the electron case where the radiative 
corrections are large it is perhaps, however, of interest to investigate the effect of 0(a 2 ) 
and higher order corrections. 

Independently of any specific conjectures on possible sources of systematic effects, 
the systematic error of 0.0045 (3.3%) assigned to the LEP average value of A T in Ref.[l] 
would appear to be unduly optimistic, in view of the apparent inconsistencies between 
measurements from both different decay channels and different experiments, discussed 
above. The value of A T derived from the euu channel deviates from the overall weighted 
average by 8.6 times the above systematic error estimate, and the L3 experiment by 4.4 
times. The Student's t common mean test shows that the uncorrelated systematic errors 
on the individual experiments are probably, on average, overestimated. Thus the real 
inconsistencies between experiments and decay channels are probably larger than those 
estimated from the assigned experimental errors. This argues even more strongly that 
the assigned systematic error on the average value of A T is too small. 

To illustrate the effect of an under-estimation of the systematic error in the r-polarisation 
measurements, consider the effect of doubling the systematic error on the LEP aver- 
age value of A T , and assigning the same systematic error to A e measured from the r- 
polarisation asymmetry^. 

This results in a weighted average value of Ai of 0.1518(26), which lies almost mid- way 
between the 'r-polarisation -in' and 'r-polarisation -out' values shown in Table 2. 

The conclusion of this discussion of the r-polarisation measurements is that there is 
evidence that the uncorrelated systematic errors of the different experiments for the same 
decay channels have been overestimated. On the other hand correlated systematic errors 
between different decay channels seem to have been underestimated. The true world- 
average value of Ai probably lies between the two values quoted in Table 2. 



The SLD A LR Measurement 

The discussion of this measurement can be very brief, as there is only one experiment 
and the present error is strongly statistics-dominated p3 |. The relative systematic error 



on A° LR = A e , which is almost completely determined by that on the SLC electron beam 
polarisation, is only 0.7%, as compared to the relative statistical error of 2.8%. The 

15 W.Lohmann, private communication 

16 In Ref.[l] a systematic error of only 0.002 is assigned to A e . Why this should be less than half of 
the already small error assigned to A T is not clear to this writer, taking into account the close similarity 
of the two measurements. Table 20 of Ref.[l] contains, however, the disclaimer that: 'the determination 
of the systematic part of each error is approximate'. 
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accuracy of the beam polarisation measurement has been checked at the 4% level by 
independent determinations using Compton and M0ller polarimeters. For comparison, 
the difference between the LEP r-polarisation and the SLD Am measurements of Ai 
amounts to 8%. The systematic error of the Am measurement would have to be wrong by 
an order of magnitude to account for this difference. Another very important advantage of 
the Alr measurement is the low level of the QED radiative correction, of only ~ —3 x 10~ 4 
24|j as compared with -110 % for A°/ B or -8% to +20%, depending on the decay channel 



and the cuts on the final state radiation |[20|| , for (P T ). 



To summarise the above detailed discussion of the different measurements contribut- 
ing to the world average value of A\ used in the model independent analysis of Section 
2; although the two LEP measurements {A°/ B and r-polarisation ) and the SLD Am 
measurement have, currently, similar statistical errors the situation is very different with 
respect to systematic errors. Only for the Am measurement is the systematic error negligi- 
ble. For A°p B and (P T ) the estimated statistical and systematic errors are almost equal |TJ. 
However, inconsistencies at the 2-3 standard deviation level in r-related measurements 
of both A°f B and (P T ) indicate that the true systematic error is probably considerably 
larger. 



Physical Interpretation of the Measured Effective 
Weak Coupling Constants 



The test of the SM provided by measurements of Z decays at LEP and SLD is, es- 
sentially, that of the SM predictions for the quantum corrections, arising from massive 
virtual particle loops, to the Born level diagrams for e + e _ — > Z — > ff. These corrections 
may be conveniently expressed in terms of two parameters Apf and Akj for each fermion 
flavour p5| . The parameters are given, in terms of the effective couplings, by the relations: 



Ap f = -2(1-2|5/|) (4.1) 
A* = ^=5^-1 (4.2) 

Ak c = 3 (1 ~ rc) - 1 ( 4 - 3 ) 



8..s 



2 
W 



4s w 

Here, following the usual on-shell definition f26|: 



AK b = 3 ( \ 2 rb) -1 (4.4) 



M 2 

sm 2 9 w = s\ v = \-c% / = \--^ [ (4.5) 

Since AKf is determined by Tf, only the weak theoretical assumption of lepton universality 
is needed to extract it from the experimental measurements. 

The SM predictions of Section 2 used the fixed values: m t = 172 GeV, m H = 149 GeV 
found in the global fit of Ref. [1]. The effect on the SM prediction of varying m t and m H 
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within the existing experimental bounds [27|, p8[ is now considered. The dependence of 



Apf on m t and ran is contained in the terms |p5 



8^ 



7T Z 



Ap f yy = ^ tan W 



11 / /m fl 



5 
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(4.6) 
(4.7) 



where £/ = for / ^ b and -4/3 for / = b. The quantum correction Akj is calculated 
using a parameterisationj^] of the ZFITTER [^{| prediction of the effective leptonic weak 
mixing angle: 



0.233597-8.95 x 10~ 8 m f 2 - 3.86 x 10~ 4 lnm t + 5.43 x Hr 4 lnm H (4.8) 



where m t and m# are in GeV units. Akj is related to 



(si 



rj)/4 by Eqns. 



(2.1), (2. 9), (2. 10) and (4.2-4.4). For the b quark there is an additional non-universal 
contribution: 



4\/27r 2 ' 



(4.9) 



The values of Apy, A/tj for / = I, c, b, extracted from the measured effective couplings 
using Eqns.(4.1)-(4.5), are presented in Table 14. Standard Model predictions are shown 
for the cases m t = 172 GeV, m H = 149 GeV and for m t = 180 GeV, m H = 100 GeV. 
The latter choice gives a somewhat better description of the leptonic corrections. In 
Table 15 corresponding results for Apf and Akj are shown for the case when the r- 
polarisation measurements are excluded from the LEP average value of A\. 

Good agreement with the SM is seen for leptons and c quarks. For b quarks however, 
the measured values of the quantum corrections are much larger than the SM predictions. 
For Apb the measured value exceeds the SM prediction by a factor of 13-15, and is of 
opposite sign. The measured value of Anb has the same sign as the SM prediction, but is 
9-11 times larger. Both effects are at the > 3 standard deviation level, but they are highly 
correlated. The discrepancies seen are so large that the significance of the deviations shows 
almost no sensitivity to m t and r/ig. 

It is also instructive to present the quantum corrections in terms of the 'epsilon pa- 
rameters' introduced by Altarelli et al. [[Jl], [H], |3"2]| . In terms of the variables used in the 
present paper to describe the effective couplings, these are defined as PHI : 



ei = Api = - 
e 2 = c 2 a Api + 



2(1 + 2a,) 
s 2 Ar w 



{cl - si) 
e 3 = c 2 A Pl + (c 2 - s 2 )Ak' 

here Sq = 1 — cl and Ar w are defined by the relations: 



2s 2 Q Ak' 



(4.10) 
(4.11) 
(4.12) 



1 



Ml. 



Ml 



2 „2 



SnC, 



0MJ 



na(M z 



1 - Ar w x/^Mf (1 - Ar w ) 



17 The relative accuracy of the formula (4.8) is about one per mille for the interesting range of values 
of mt and m# . 
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and 



In Ref. [32] a fourth parameter, e&, was introduced. It may be defined in three distinct 
ways: 

e b (a b ) = ^-1 (4.13) 
a; 

= (4.14) 

l - n 

e ' (Sl) = 2(*)»(1-6 W + 2K,) (4 ' 15) 
3 

In the SM, retaining only the leading terms ~ mf, the three definitions (4. 13)- (4. 15) are 
equivalent! 1 ^. In previous phenomenological applications however, |32|, Q only the third 



where 



definition (4.15) based, via Eqn.(2.8) on the measurement of R b was used. The measured 
values of the six epsilon parameters defined above are presented in Table 16, where they 
are compared with the SM predictions. As noted previously |jvj the values of ei, e 2 and e 3 
are in good agreement with the SM predictions. A small deviation is observed for e b (s b ), 
a residual of the much discussed |l[ l R b problem'. However, both e b (a b ) and e b (f b ) deviate 
from the SM prediction by about four standard deviations^]. One may note the extreme 
sensitivity of the parameter e b (r b ) to the anomalous b coupling; the measured value is 39 
times and 4.7a0 larger than the SM prediction. The SM predictions for these quantities 
are insensitive to vtlh and are essentially given by the term ~ m 2 : 



2 Cm 2 

e b = --Ap top = = -0.0062 (m t = 172GeV) (4.16) 

3 4v27T 2 

Table 17 shows the epsilon parameters calculated excluding the r-polarisation measure- 
ments from the averages. It can be seen that the measured value of e b (f b ) exceeds the SM 
prediction by a factor of 44 and more than six standard deviations in this case. It may be 
remarked that, in this case, the leptonic parameter e 3 also shows a deviation of 3 to 4cr. 

The conclusion to be drawn from Tables 16 and 17 is that the deviations observed for 
the b quark couplings, interpreted as a real physical effect, do not enter at all into the 
framework of the SM nor any of its 'natural' extensions. Supersymmetry, Technicolour, 
anomalous WW^ or WWZ couplings, and new U(l) gauge bosons are all expected, via 
vacuum polarisation effects in the gauge boson propagators, to produce deviations from 
the SM predictions for ej, e 2 or e 3 pp| , |3~T| , |3"2] ]. These parameters are much more sensitive 



to any flavour independent modifications of the couplings, due to anomalous vacuum 
polarisation effects, because of the high precision of the purely leptonic measurements. 



18 Modulo small b-mass dependent corrections 

19 Again, the errors on these quantities are highly correlated 

20 The errors on this quantity, determined essentially by those on A b , are skewed and non-gaussian. 
The average error is quoted in Tables 16 and 17. The confidence level of the deviation of e b (r b ) from the 
SM, assuming gaussian errors for A b , is in fact almost the same as that of the latter, about one per mille. 
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Although the €3 parameter shows a quite large deviation from the SM prediction in the 
case that the r-polarisation measurements are excluded, the most important apparently 
anomalous effect occurs in the quantum corrections to the b quark couplings that disagree, 
by an order of magnitude, with the expectations of the SM. 

A clue as to the origin of the anomalous b quark couplings is provided by considering 
the right- and left-handed effective couplings, g^, g% related to a b and v b by the relations: 

9b=l(v b -a b ) = -^p~ b e b {s b w ) 2 (4.17) 
9 b =\(v b + a b ) = ^T b [T b 3 - e b (s h w f] (4.18) 

From the measured values of a b and v b presented in Table 7, the following values of the 
left-handed and right-handed effective couplings of the b quarks are found: 

g\ = -0.4155(30) #f = 0.1098(101) 

which may be compared with the SM predictions of: 

g\ = -0.4208 #f = 0.0774 

The value of g b is quite consistent with the SM prediction ( for m t = 172 GeV, m B = 149 
GeV ) ( a 1.3%, 1.8a deviation) whereas the discrepancy for g b is much larger, (a 42%, 
3.2a deviation). Excluding the r-polarisation measurements gives the results: 

#£ = -0.4138(29) gf = 0.1160(90) 

The deviations from the SM prediction are 1.7% and 2.4a for gfc and 50% and 4.3a for 
~g~b '. One may remark that the weak iso-spin of the SM affects only g£ not g^, so it is 
possible that the SM does correctly describe g% but that there is a new, anomalous, right 
handed coupling for the b quark. 

The right- and left-handed effective couplings of the s,d quarks have recently been 
measured by the OPAL collaboration [Q with the results: 



jjL _ _n 44 +0.i3 -R _ n 1 0+0.15 
9d,s — u - 44 -0.09 9d,s — u - l6 -0.l7 

to be compared with the SM predictions -0.424 and 0.077 respectively. These measure- 
ments are in good agreement with both the SM predictions and the measured b quark 
couplings given above. 

Limits can also be set on possible anomalous couplings of the other 'd-type' quarks, 
d,s by comparing the measured values of (A 9 FB ) and Thad with the predictions of a 
model in which the d and s quarks are assumed to have the same effective coupling 
constants as those measured for the b quarks. The prediction of this model for (A q FB ) is 
0.1600(72), which is consistent with the 'measured' value (see Section 2 above and Table 
8) of 0.1592(86) at the 0.68a level. No useful constraint on possible anomalous couplings 
of the d and s quarks, of a size similar to those observed for the b quark, is therefore 
obtained using the measured value of (A q FB ) with the present experimental errors. A 
more favourable case is T had . Using the world average value of a s (M z ) of 0.118(5) || [U| 
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in Eqn.(2.17) to calculate the QCD correction, and with C® ED = 1.00040, the predicted 
value of Thad in the model with a universal right-handed anomaly for down-type quarks is 
1.7249(46) GeV. This differs by 3.6<r from the LEP average measurement |T| of 1.7436(25) 
GeV. Thus the model is essentially excluded by the measurement of T had . It is interesting 
to note that the precise measurement of currently gives a much more stringent con- 
straint on possible anomalous couplings of the d and s quarks than the direct measurement 
of their left- and right-handed couplings cited above |34] . 



The values of the left- and right-handed couplings of the c quarks derived from the 
measured values of a c and v c given in the first two rows of Table 7 are: 



g L c = 0.3440(92) g? = -0.1600(70) 

in very good agreement with the SM predictions for m t = 172 GeV, m# = 149 GeV of: 

g L c = 0.3465 gf = -0.1545. 

The ±2<t limits for deviations of gf from the SM prediction extends from -0.174 to -0.146. 
Thus, at 95% CL, any anomalous right-handed couplings of the c quark lie between -9% 
and +15% of the SM prediction. 



5 Summary and Outlook 



The analysis presented in this paper has been performed in such a way as to separate, 
as far as possible, the actual results of experimental measurements, expressed in terms 
of the effective weak coupling constants of the leptons, c quarks and b quarks, from the 
comparison of these results with the SM. The quantities vi, ai and Tq = vq/clq ( Q =c,b ) 
were extracted assuming only lepton universality. A further theoretical assumption is 
needed to extract vq and Hq separately. The hypothesis of non-b quark lepton universality 
was chosen. It was shown that the alternative of assuming a s (Mz) to be knownQ, would 
give essentially the same results for the b quark couplings, if the world average value of 
a s (M z ) = 0.118(5) H [K| is used. Thus the actual non-b quark couplings are consistent 
with the universality hypothesis. 

The measured effective couplings of the leptons and c quarks are in good agreement 
with the SM predictions for a values of m t and m# well consistent with existing experimen- 
tal limits. However, the b quark couplings deviate from the SM predictions by more than 
three standard deviations. Taking carefully into account all relevant error correlations, 
the probablity that all six effective couplings are consistent with lepton universality and 
the SM is found to be 0.9%. This number is the product of the CL for consistency with 
lepton universality, given by the x 2 of the different Ai when compared with their weighted 
average (8.4 %) and the CL of the SM comparison using the average value from Table 
9 (10.5 %). The probability that the leptonic and b quark couplings are consistent with 
lepton universality and the SM is similarly calculated to be 0.18%. The latter probability 
drops to only 0.018% if the r-polarisation measurements are excluded in calculating the 
average value of A\. 

21 This was the procedure adopted in Ref. [7], where where heavy quark effective couplings consistent 
with those extracted in the present paper were found. 
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In Section 3 the possibility is discussed that systematic effects, beyond those taken 
into account in the present analyses and correlated between different experiments, exist 
in the r-polarisation measurements 

In Section 4 the measured effective couplings are analysed in terms of the one-loop 
quantum correction parameters Apf and AtZf. These quantities are found to be in agree- 
ment with the SM predictions for the leptons and c quarks, but to be an order of magni- 
tude larger than these predictions for the b quarks. An analysis in terms of the 'epsilon' 



parameters |30], [31], [32[ ei, 62, 63, efc(s&), eb(a&) and e&(r&) shows reasonable agreement 
with the SM for the first four parameters, but ~ 4 — 6a deviations for the last two. The 
agreement found between experiment and the SM for e±, e% and €3 demonstrates that 
there is no experimental evidence for contributions from natural extensions of the SM 
such as Supersymmetry, Technicolour, anomalous WW7 or WWZ couplings, or new U(l) 
gauge bosons, which are all expected to produce deviations from the SM model predic- 
tions for one or more of these parameters pi]] . It may be remarked, however, that if the 
r-polarisation measurements are excluded the €3 parameter differs from the SM predic- 
tion by three standard deviations, even for a Higgs boson mass (run = 100 GeV) close 



to the current experimental lower limit [28] of about 70 GeV. This additional potential 



problem for the SM merits further investigation, but is beyond the scope of the present 
paper. 

A simple picture emerges in terms of the right-handed and left-handed couplings of 
the b quark. While the latter agrees at the 2a level or better with the SM prediction, 
the former shows a 42-50% and 3.2-4.3<r deviation. The larger deviations occur when 
r-polarisation measurements are excluded. Thus the most significant deviation from 
the SM predictions of the effective couplings is found in the right-handed b quark weak 
coupling constant. 

The large deviations from the SM found here for the b quark couplings confirm results 
presented in a previous paper |7|. The paper contains, however, neither any discussion 
of the overall statistical significance of the deviations nor any physical interpretation of 
them. Other recent reviews |33], |35|, Q did not extract the heavy quark couplings from 



the measurements. They discussed rather the high precision purely leptonic data and 
the quantities R c and Rb all of which agree well with SM predictions. Of the three 
parameters €b{~5b), eb(a&) and ef,(r&) defined in Ref.[32], only e&(s&), derived from Rb, was 
considered in Ref.[33]. As also only the ALEPH Rb measurement |37|] (which is in almost 



perfect agreement with the SM) was used, the anomalies in the b quark couplings were 
undetected. As shown in Tables 16 and 17 it is the other two 6b parameters that are most 
sensitive to the apparent deviations from the SM. 

The most important remaining question is whether the observed deviation in the 
right-handed b quark coupling is likely to be confirmed or excluded by currently planned 
measurements of electroweak observables. The effect observed in the effective couplings, 
interpreted as a statistical fluctuation, is ~ 3-4cr. A similar sized deviation for Rb was 
reported in 1995 by the Electroweak Working Group [J38| . A year later, as a result of a 
better systematic understanding of both the Rb and the correlated R c measurements, as 
well as much improved statistical errors, the deviation of the average Rb measurement 
was reduced to below two standard deviations ffl. Can the presently observed anomaly 
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in the b quark right-handed coupling be expected to meet the same fate as that of 
Arguments will now be given that this is unlikely, given the current status of the analysis 
of the LEP+SLD electroweak data. If there is some large, as yet unknown, systematic 
effect in the LEP A°p b B and/or the SLD Ab measurements similar to those uncovered during 
1996 in the Rb and R c measurements, the anomaly might disappear (or become larger). 
This hypothetical question will not be further addressed here. 

The discussion concentrates on the observed (and possible future) deviations from 
the SM of Ft, which is extracted from measurements using only the weak theoretical 
assumption of lepton universality. As can be seen in Figs, lc, 2c this quantity shows much 
larger deviations from the SM predictions than s&. The measurements of several of the 
relevant electroweak observables: A°/ B , Alr, A°p b B and Ab have recently been improved as 



compared to numbers quoted in Table 1 fl39[. The updated value of following the same 



analysis procedure as described in Section 2 above is 0.591(30). This is very consistent 
with the value reported in Table 3. The discrepancy with the SM is slightly reduced 
from 3.34 to 3.26 standard deviations. The forseen improvements in the precision of 
measurements of electroweak observables have been discussed in detail in the Review of 
Renton |7J]. Only very modest improvements are to be expected from completing the 
analysis of the existing and final LEP1 data. Their impact on the Tb measurement is 
expected to be essentially negligible. As shown in Sections 2 and 3 above this is not 
true, however of possible systematic effects, witness the sensitivity of the CLs in Table 9 
and the Ap, An and epsilon parameters in Tablesl4-17,to the inclusion or exclusion of 
the r-polarisation data. A better systematic understanding of the latter could have a 
large effect on the LEP average value of Ai. Ultimately however the LEP+SLD average 
value of A[ is expected to be dominated by the improved precision of the ongoing A LR 
measurement. According to Ref. [7] the errors on the existing SLD measurements of Alr 
and Ab [|39| should be reduced, by the end of the SLD experimental program, by 41% and 
45% respectively. The most significant improvement is expected to be in the accuracy 
of the LEP+SLD average value of Ai due to the Alr measurement. The error of the 
latest measurement: Af LD = 0.1547(32) |39[ is expected to be reduced to 0.0019. This is 
significantly smaller than the error, 0.0033, on the LEP average value of Ai which will, in 
the absence of new, presently unknown, systematic corrections, change little in the future. 
The improvement on the most recent SLD measurement of Ab, A^ LD = 0.897(47) ||39|| , 
where the error is expected to be reduced to 0.0027, is also large. However, the LEP 
average value of A b , derived from A FB is already more precise: A LEP = 0.861(21). To 
study the possible impact of the new SLD measurements on the expected future value of 
Tb the value of Af LD is allowed to vary by ±2a G (where o"o = 0.0019 is the final expected 
error) about the latest measured value. The weighted average of Af LD and the latest 
LEP measurement, is calculated (including, or not, the r-polarisation data) to estimate 
the likely range of the 'final' LEP+SLD average value of A[. For each Ai value, Ab is 
extracted from the LEP average A°p b B using Eqn.(2.3) to yield 'final' values of A^ EP . The 
weighted average of A LEP and A^ LD is then made, assuming for the latter the latest value 
given above and the expected final error of 0.0027. In the last step, for each 'final' Ab 
value Tb is extracted using Eqn.(2.6) and compared with the SM prediction. The results 
of this exercise are presented in Table 18. This purely statistical study shows that the 
deviations in Tb are unlikely to drop below 3a , or to become more significant than 5a. 
Larger deviations are still seen when the r-polarisation data are excluded, but the effect 
is less than shown in Table 3 (~ 0.5<r instead of ~ la) due to the much higher statistical 
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weight expected from the final SLD A\ measurement. Unlike in the case of the Rt, anomaly 
in 1995, most of the foreseeable LEP+SLD data contributing to the measurement of Fj, 
have already been analysed. The predictions presented in Table 18 show that only small 
changes are to be expected in the statistical significance of the observed anomaly in the 
b quark couplings. 

The most important message of this paper for future electroweak analyses is that the 
step of extracting the important physical quantities from the experimental measurements 
should be clearly separated from the comparison of these quantities with theoretical pre- 
dictions. This was done routinely in the past for the leptonic weak coupling constants |]]. 
A similar procedure should be followed, in future, also for the heavy quark couplings. 
Only in this way can the physical origins of apparent deviations from theory be precisely 
located, and the CLs for agreement of the measurements with the theory be easily and 
correctly calculated. The second point (exemplified here by the discussion of the LEP 
r-polarisation measurements) is to check carefully the internal consistency of different 
experimental measurements of the same physical quantity before they are averaged. This 
can give indications of hitherto unconsidered systematic effects. In any case, the CL for 
the internal consistency of different measurements used to calculate an average physical 
quantity should be multiplied by that given by the x 2 of the theory/experiment com- 
parison based on the averaged quantities, in order to give a more meaningful overall 
probability for the theory/experiment comparison. 
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Table 1: Average values of electroweak observables used in the analysis IQ. SM denotes 
the Standard Model prediction for m t = 172 GeV, m# = 149 GeV |1]] 



A l 


A c 


A b 


R c 


Rb 


0.1501(24) 
0.1533(27) 


0.645(39) 
0.634(38) 


0.869(22) 
0.853(22) 


0.1715(56) 


0.2176(11) 



Table 2: LEP+SLD averages. The r-polarisation measurements are excluded from the 
averages quoted in the second row 







r c 


r b 


Measurement 
SM 

(Meas.-SM) /Error 


0.07548(120) 
0.07332 
-1.80 


0.366(29) 
0.383 
-0.59 


0.582(32) 
0.689 
-3.34 


Meas. r poln. out 
(Meas.-SM) /Error 


0.07711(140) 
-2.71 


0.357(29) 
-0.90 


0.562(29) 
-4.38 



Table 3: Measured values of r f = Vf/cif compared to Standard Model predictions 





Sl 


S c 


s b 


Measurement 


0.25244(33) 


0.2877(95) 


0.3676(24) 


SM 


0.25259 


0.2880 


0.3644 


(Meas.-SM) /Error 


-0.45 


-0.03 


1.33 



Table 4: Measured values of Sf = aj(l — 6fif) + v compared to Standard Model 
predictions 
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leptons 




Meas. 


SM 


Dev(o-) 




-0.50101(33) 


-0.50124 


0.67 


Vl 


-0.03782(68) 


-0.03675 


-1.57 


ai 


-0.50093(33) 




0.91 


Vl 


-0.03863(77) 




-2.44 



Table 5: Measured values of the effective electroweak coupling constants for the charged 
leptons. Dev(a) = (Meas.-SM )/Error. The values given in the last two rows exclude 
r-polarisation data from the averages 



C QED 




(JQCD 




1.00046 


0.99975 


1.0012 


0.9953 



Table 6: QED and QCD correction factors for heavy quarks assuming a s (Mz) = 0.12 
and a(M z )~ l = 128.9. 





c 


quark 




b quark 




Meas. 


SM 


Dev(o-) 


Meas. 


SM 


Dev(o-) 


a f 


0.504(10) 


0.501 


0.30 


-0.5252(75) 


-0.4981 


-3.61 




0.184(15) 


0.192 


-0.53 


-0.3057(125) 


-0.3434 


3.18 


a f 


0.505(10) 




0.40 


-0.5298(70) 




-4.53 




0.180(15) 




-0.80 


-0.2977(123) 




3.72 



Table 7: Measured values of the effective electroweak coupling constants of c and b 
quarks. Dev(a) = (Meas.-SM)/Error. The values given in the last two rows exclude 
r-polarisation data from the averages 





'Measured' 


SM Pred. 


MI Pred. 


MI Pred. with SM b quark 


(A q FB ) 


0.1592(86) 


0.1641 


0.1639(28) 


0.1692(28) 


('Meas'.-Pred.)/Error 




-0.57 


-0.52 


-1.1 



Table 8: Values of the mean quark charge asymmetry. 'MI Pred.' stands for Model 
Independent Prediction (see text). See also the text for the definition of 'measured'. 



Observables 


Th Si 


ri, n 


n, si, r b , s b 


ri, r b , r c 


r h si, r b , s b , r c , s c 


dof 

x 2 

CL (%) 


2 
3.44 
17.9 


3 
10.6 
1.4 


5 
13.2 
2.2 


5 
10.9 
5.3 


8 
13.2 
10.5 


x 2 

CL (%) 


7.55 
2.3 


17.2 
0.064 


19.4 
0.16 


17.4 
0.38 


19.6 
1.2 



Table 9: x 2 an d confidence levels for agreement with the SM ( m t = 172 GeV, uih = 
149 GeV) of different sets of electroweak observables sensitive to the effective couplings, 
assuming perfect statistical consistency of the LEP+SLD averages in Table 2. The values 
given in the last two rows do not use r-polarisation measurements in the A\ average. See 
the text for the explanation of the number of degrees of freedom (dof) in each case. 
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T — fJ, 


t — e 


fj, — e 


CL(%) STT 


2.7 


6.0 


68 


CL(%) EET 


8.0 


17 


94 



Table 10: Confidence levels for the consistency of LEP measurements of A FB , Ap B and 
A T p° B . STT: Student's t test, EET: CL calculated from Estimated Experimental errors. 





ALEPH 


DELPHI 


L3 


OPAL 


WA 


evv 
pvv 

7TV 

pv 
a\v 


0.200(51,31) 
0.124(41,21) 
0.142(20,11) 
0.108(19,18) 
0.135(35,20) 


0.179(52,67) 
0.097(38,22) 
0.158(33,50) 
0.199(46,39) 
0.103(50,38) 


0.168(39,15) 
0.111(45,16) 
0.135(21,17) 
0.168(17,10) 


0.161(33,29) 
0.138(33,22) 
0.117(14,12) 
0.116(13,11) 
0.151(37,31) 


0.173(26) [16] 
0.119(22) [11] 
0.130(12) [9] 
0.135(12) [22] 
0.134(28) [14] 


WA 


0.132(15) [24] 


0.137(26) [27] 


0.154(14) [14] 


0.123(11) [9] 


0.1343(73) 



Table 11: LEP measurements of A T = —(P T ) The first error shown is statistical, the 
second systematic. For the Weighted Averages (WA) the weighted total error is given, 
where statistical and systematic errors are added in quadrature. Errors calculated from 
sample variances are shown in square brackets. 





e — p 


e — 7r 


e-p 


e — ai 


/! — 7T 


p- p 


p — a\ 


7T — p 


TX — OL\ 


p-o-i 


CL(%) STT 


0.3 


1.5 


36 


2.5 


15 


35 


9 


68 


45 


95 


CL(%) EET 


11 


13 


19 


31 


67 


51 


66 


74 


79 


97 



Table 12: Confidence levels for the consistency of LEP (P T ) measurements using different 
r decay channels. STT and EET are defined in Table 10. 





A-D 


A-L 


A-0 


D-L 


D-O 


L-0 


CL(%) STT 


85 


31 


15 


31 


31 


64 


CL(%) EET 


87 


29 


63 


57 


63 


9 



Table 13: Confidence levels for the consistency of LEP (P T ) measurements by different 
experiments. STT and EET are defined in Table 10. 





leptons 


c quark 


bq 


uark 




Ak, 


Ap c 


Ak c 


Ap 6 




Expt. 


0.00404(133) 


0.03445(134) 


0.016(41) 


0.064(49) 


0.101(32) 


0.403(107) 


SM m t = 172 GeV 
m H = 149 GeV 
Dev(cr) 


0.00497 
-0.7 


0.03686 
-1.8 


0.005 
0.27 


0.037 
0.55 


-0.007 

3.38 


0.0436 
3.36 


SM m t = 180 GeV 
m H = 100 GeV 
Dev(a) 


0.00563 
-1.2 


0.03472 
-0.02 


0.006 
0.24 


0.034 
0.61 


-0.008 
3.40 


0.0412 
3.38 



Table 14: Measured values of the quantum correction parameters Apf and A/tj compared 
to SM predictions. Dev(c) = (Expt.-SM)/Error. 
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lept 


ons 


c quark 


b quark 








Ap c 


Ak c 


Ap 6 


An b 


Expt. 


0.00372(133) 


0.03260(157) 


0.020(41) 


0.079(49) 


0.119(30) 


0.470(97) 


SM m t = 172 GeV 














m H = 149 GeV 


0.00497 


0.03686 


0.005 


0.037 


-0.008 


0.0439 


Dev(cx) 


-0.94 


-2.7 


0.37 


0.86 


4.23 


4.39 


SM m t = 180 GeV 














m H = 100 GeV 


0.00563 


0.03472 


0.006 


0.034 


-0.008 


0.0412 


Dev(a) 


-1.4 


-1.4 


0.34 


0.92 


4.23 


4.42 



Table 15: Measured values of the quantum correction parameters Apf and Akj com- 
pared to SM predictions, r-polarisation measurements are excluded from the A t average. 
Dev(cr) = (Expt.-SM)/Error. 





ei 


£2 


£3 




e&K) 




Expt. 


0.00404(133) 


-0.0073(8) 


0.0031(8) 


-0.0017(18) 


0.048(15) 


-0.263(57) 


SM m t = 172 GeV 














m H = 149 GeV 


0.00497 


-0.0076 


0.0051 


-0.0045 


-0.0060 


-0.0068 


Dev(cr) 


-0.7 


0.38 


-2.5 


1.6 


3.6 


-4.5 


SM m t = 180 GeV 














m H = 100 GeV 


0.00563 


-0.0062 


0.0045 


-0.0046 


-0.0068 


-0.0055 


Dev(a) 


-1.2 


-1.4 


-1.8 


1.6 


3.7 


-4.5 



Table 16: Measured values of the epsilon parameters of Refs. [16-18] compared to SM 
predictions. Dev(a) = (Expt.-SM)/Error. 





ei 


£2 


£3 


£b(s b ) 


efe(ab) 




Expt. 


0.00372(133) 


-0.0068(8) 


0.0021(8) 


-0.0018(18) 


0.058(14) 


-0.298(47) 


SM m t = 172 GeV 














m H = 149 GeV 


0.00497 


-0.0076 


0.0051 


-0.0045 


-0.0060 


-0.0068 


Dev(cr) 


-0.94 


1.0 


-3.8 


1.5 


4.6 


-6.2 


SM m t = 180 GeV 














m H = 100 GeV 


0.00563 


-0.0062 


0.0045 


-0.0046 


-0.0068 


-0.0055 


Dev(cx) 


-1.4 


-0.75 


-3.0 


1.6 


4.6 


-6.2 



Table 17: Measured values of the epsilon parameters of Refs. [30-31] compared to SM 
predictions, r-polarisation measurements are excluded from the Ai average. Dev(o") = 
(Expt.-SM)/Error. 
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'final values' of Ai 


'final values' of r b 


A? LD ±x,x = 


-2<7 





+2a 


-2cr 





+2a 


LEP+SLD 
Full average 


0.1501(16) 


0.1525(16) 


0.1553(16) 


0.604(26) 
[-3.3] 


0.592(25) 
[-3.9] 


0.579(23) 
[-4.8] 


t poln. 
Excluded 


0.1513(17) 


0.1545(17) 


0.1577(17) 


0.598(26) 
[-3.5] 


0.583(24) 
[-4.4] ' 


0.568(23) 
[-5.3] ' 



Table 18: A scenario for the 'final' values and errors of Ai and r b at the end of the 
LEP+SLD experimental programme. Af LD = 0.1547(19) and Af LD = 0.897(27) are 
assumed, and Af LD varies by -2a to +2a from the central value, where a is its error. 
For f b) the number in square brackets denotes: (Value-SM)/Error, where rf M = 0.689. 
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